A third massive star component in the a Orionis AB system 
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ABSTRACT 

We report on the detection of a third massive star component in the a Orionis AB system, 
traditionally considered as a binary system. The system has been monitored by the IACOB 
spectroscopic survey of Northern Massive stars program, obtaining 23 high-resolution FIES@NOT 
spectra with a time-span of ~2.5 years. The analysis of the radial velocity curves of the two 
spectroscopic components observed in the spectra has allowed us obtain the orbital parameters 
of the system, resulting in a high eccentric orbit (e ~ 0.78) with an orbital period of 143.5 ± 0.5 d. 
This result implies the actual presence of three stars in the a Orionis AB system when combined 
with previous results obtained from the study of the astrometric orbit (with an estimated period 
of ~ 157yr). 



Subject headings: stars: early-type — binaries: spectroscopic — stars: individual (a Ori AB) 
star clusters: individual (a Orionis) 



- galaxies: 



1. Introduction 

The a Ori star (V ~ 3.80 mag, HD 37468, 
HR 1931) is visible with naked eye at about 1 deg 
to the south of Alnitak (£ Ori), the easternmost 
bright star of the Orion Belt, and at a compara- 
ble separation to the west of the Horsehead Neb- 
ula. The actual multiple status of a Ori was 
firstly described in Tabula Nova Stellarum Du- 
plicium, where Mayer (1779) tabulated the stars 



1 Based on observations made with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 



known today as a Ori A and a Ori E, separated 
by 42 arcsec, and proposed a third component in 
between. Struve et al. (1876), from observations 
between 1819 and 1831, confirmed the third com- 
ponent, a Ori D, and reported a new "ash" -colored 
one, a Ori C. The brightest star in the system 
was found to be, in its turn, a tight binary by 
Burnham (1892). The two components, histori- 
cally known as a Ori A and a Ori B, are thought to 
have very early spectral types (09.5 V and BO. 5 V, 
respectively) and are separated by no more than 
0.3 arcsec. 

The angular separations p and 6 between 
a Ori A and B have been frequently monitored 
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during the last 120 years by numerous authors 
with micrometers, first, and speckle imaging and 
adaptive optics, next. Some orbit determinations 
have been published by Kiimmritz (1958), Heintz 
(1974, 1997), Hartkopf et al. (1996), and Turner 
et al (2008). The relatively long orbital period 
(P ~ 157 yr; i.e., the binary has not completed 
a whole revolution since its discovery), wide pro- 
jected physical separation (s ~ 100 AU, assuming 
d w 350-450 pc) and early spectral types have sug- 
gested that a Ori A and B may form "the most 
massive visual binary known", with up to 40 Mq 
(Heintz 1974; Caballero 2008a). 

The small orbital semimajor axis of only 
0.266±0.002 arcsec and nearly circular orbit (Turner 
et al. 2008) makes impractical an individu- 
alized spectroscopic study. Predicted radial- 
velocity variations in the primary have a consid- 
erable peak-to-peak amplitude of about 7kms _1 
(Hartkopf et al. 1996), but the long orbital pe- 
riod, magnitude difference between both compo- 
nents (AHp = 1.21±0.05mag; Perryman et al. 
1997) and scarcity of adequate high-resolution 
spectra in the last century complicate any long- 
term spectroscopic study. Some authors have re- 
ported, however, that at certain epochs the spec- 
trum of a Ori "A— B" changes in appearence from 
its single-line stage, showing broadened and dou- 
ble lines separated by up to 280 km s -1 (Frost & 
Adams 1904; Frost et al. 1926; Miczaika 1950; 
Bolton 1974; Fullerton 1990; Morrell & Levato 
1991; D. M. Peterson, priv. comm.). The time 
scale of these variations was only a few months. 
Both observational results would not be expected 
from the study of the astrometric orbit, and could 
suggest that the a Ori "A-B" system is a triple 
system (one of the stars resolved visually being 
a double-lined spectroscopic binary, viz. Bolton 
1974). Nevertheless, the vast majority of the spec- 
trophotometric observations of a Ori "A-B" have 
not been able to resolve the spectroscopic binary 
and, in words of Frost & Adams (1904), "evi- 
dences of complexity in the spectrum [of a Ori 
are] scarcely sufficient to justify conclusions on 
the subject". More than one century later the 
"subject" is not resolved yet. 

The Trapezium-like a Ori system being a dou- 
ble or a triple massive star has major implications 
because of its location in the center of the homony- 
mous a Orionis cluster, which is a cornerstone for 



star formation studies (Garrison 1967; Wolk 1996; 
Bejar et al. 1999; Walter et al. 2008; Caballero 
2008b). First, the top of the (initial) mass func- 
tion and the long-term dynamical evolution of the 
cluster depends on the actual mass of the bright- 
est stars (Sherry et al. 2004; Caballero 2007). 
Second, any derivation of theoretical masses of 
intermediate- and low-mass stars and substcllar 
objects rely on a precise value of the cluster helio- 
centric distance, which is often determined from 
fits of the cluster sequence in color-magnitude di- 
agrams (e.g., Sherry et al. 2008; Mayne & Naylor 
2008) or via dynamical parallax (Caballero 2008a); 
both procedures must take into account the multi- 
plicity of a Ori "A-B" . Finally, the "photo-erosion 
of pre-existing cores" (Whitworth et al. 2007 and 
references therein) by the strong high-energy radi- 
ation emitted by the Trapezium-like system may 
partially explain the copiousness of cluster brown 
dwarfs and planetary-mass objects (Zapatero Os- 
orio et al. 2000; Bejar et al. 2001; Caballero et al. 
2007; Lodieu et al. 2009; Bihain et al. 2009), 
not counting the nature of the nearby proplyd 
a Ori IRS1 AB (van Loon & Oliveira 2003; Ho- 
dapp et al. 2009) or the shape of the Horsehead 
Nebula (Abergel et al. 2003; Pound et al. 2003; 
Habart et al. 2005). To sum up, any investigation 
of the a Orionis cluster is incomplete if the central 
a Ori star system is not correctly accounted for. 

2. Observations, analysis and results 

The a Ori "A-B" system is one of the mas- 
sive binary/multiple systems monitored by the IA- 
COB spectroscopic survey of Northern Galactic 
OB stars project (Simon-Diaz et al. 2011). Using 
the cross-dispersed echelle spectrograph FIES at 
the Nordic Optical Telescope, we have gathered so 
far 23 high- resolution (R = 46,000) spectra with a 
time-span of ~2.5 years (between November 2008 
and April 2011; see Tabled]). Our observations 
have been complemented with one spectrum each 
of a Ori D (B2 V) and a Ori E (B2 Vpe), the other 
two massive star components in the a Ori system. 
We refer to Simon-Diaz et al. (2011) for a de- 
scription of the spectroscopic observations and the 
IACOB database. 

Representative examples of the obtained spec- 
tra are shown in Fig. [1] along with the spectra 
of a Ori D and E for comparison. At first glance, 
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Table 1: Radial velocity measurements of a Ori Aa 
and Ab. 



Date 


HJD 


^radjAa 


Urad.Ab 


(yyyy/mm/dd) 


(-2400000) 


(knis" 1 ) 


(kms" 1 ) 



2008/11/05 


54776 


736 


50.1 ± 3.5 


6.0 ± 1.9 


2008/11/05 


54776 


745 


50.2 ± 3.6 


7.0 ± 1.9 


2008/11/06 


54777 


697 


50.4 ± 3.6 


6.2 ± 2.0 


2008/11/06 


54777 


699 


50.3 ± 3.6 


6.3 ± 2.0 


2008/11/06 


54777 


701 


50.4 ± 3.5 


5.9 ± 1.9 


2008/11/07 


54778 


718 


49.1 ± 3.5 


8.0 ± 1.9 


2008/11/07 


54778 


720 


49.3 ± 3.5 


7.6 ± 2.0 


2008/11/08 


54779 


704 


49.1 ± 3.5 


6.3 ± 2.0 


2008/11/08 


54779 


706 


46.9 ± 3.6 


7.6 ± 1.9 


2009/11/09 


55145 


662 


18.2 ± 3.5 


48.0 ± 1.9 


2009/11/11 


55147 


650 


14.4 ± 3.4 


52.6 ± 1.9 


2010/09/07 


55447 


731 


-52.6 ± 1.2 


128.2 ± 1.0 


2010/09/09 


55449 


747 


-80.3 ± 1.5 


167.5 ± 1.2 


2010/10/23 


55493 


735 


45.1 ± 3.5 


0.6 ± 1.9 


2010/10/23 


55493 


737 


45.8 ± 3.4 


0.1 ± 1.9 


2011/01/11 


55573 


509 


18.5 ± 3.6 


40.8 ± 2.0 


2011/01/15 


55577 


506 


11.3 ± 3.5 


50.0 ± 1.9 


2011/01/15 


55577 


509 


12.5 ± 3.4 


50.0 ± 1.9 


2011/01/15 


55577 


512 


11.9 ± 3.5 


50.2 ± 2.0 


2011/02/11 


55604 


368 


34.4 ± 3.5 


28.3 ± 2.0 


2011/02/20 


55613 


359 


50.3 ± 3.3 


9.3 ± 2.0 


2011/03/27 


55648 


379 


52.5 ± 3.4 


5.1 ± 1.9 


2011/04/08 


55660 


353 


47.1 ± 3.4 


7.8 ± 2.0 



the spectra of a Ori "A-B" show two spectroscopic 
components, a narrow one and a broad one, vary- 
ing on a time-scale of weeks. The shown spec- 
tra have been corrected for heliocentric velocity 
and shifted to the systemic velocity of the cluster 
(+30.93 ± 0.92 kms -1 ), as determined by Sacco et 
al. (2008) from dozens of FLAMES spectra of sin- 
gle low-mass stars in the a Orionis cluster. While 
the spectra of a Ori D and E show single, broad 
components at the cluster systemic velocity, those 
of the two spectroscopic components of a Ori "A- 
B" move around that value. 

The He IA5875A line was used to determine 
the radial velocity of both spectroscopic compo- 
nents of a Ori "A-B" for the 23 compiled spec- 
tra. This strong line makes easier disentagling 
both components and provides more reliable ra- 
dial velocity measurements than other metal lines, 
for which the broad component appears too faint. 
Besides, the helium line is close to the interstellar 
Na I AA5890,5895 A doublet, which can be used as 
a sanity check of velocity shifts and heliocentric 
velocity correction. 



1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Fig. 1. — FIES spectra of a Ori Aa+Ab+B at five rep- 
resentative epochs (top) and of the B2 V stars a Ori D 
and E (bottom) around the He I A5875 Aline, corrected by 
systemic velocity and separated by 0.25 arbitrary units in 
normalized flux. The sharp, deep absorption lines are in- 
terstellar Na I AA5890,5895 A. The other shallow absorption 
lines are telluric. 



The radial velocities were determined by means 
of a two-parameter cross correlation of the ob- 
served spectra with a grid of synthetic spec- 
tra built with two rotationally-broadened, radial 
velocity-shifted (u ra d,i and u ra d,2) He IA5875A 
line-profiles computed with the fastwind stellar 
atmosphere code (Puis et al. 2005). Previously, 
the projected rotational velocities had been de- 
termined by applying the Fourier method (Gray 
1976; see also Simon-Diaz & Herrero 2007 for a 
recent application to OB stars) to the spectrum 
observed on 2010/09/09, when both components 
were not blended. We derived vsini values of 
30 ±3 and 130 ±10 km s -1 for the narrow and 
broad components, respectively. 

The resulting radial velocities, together with 
the associated uncertainties, are listed in Table Q] 
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Table 2: Orbital parameters of a Ori Aa and Ab. 



o 



i 5 
? 
- -5 

< 5 



5 



: <t(0-C)= 1-9 km s~' 


A 




1 J 


\ **** 


a a 




& : 


: a(0-C)= 2.2 km s" 1 


□ 




□ ~ 

-6--n J 











0.0 O.b 1.0 

(T = 2455452.03) 



Fig. 2. — Upper panel: Radial velocity curves of 
a Ori Aa (triangles) and Ab (squares) phased to 
the period P = 143.5 d. Errorbar sizes are smaller 
than those of symbol sizes. Vertical dotted line in- 
dicates the systemic velocity (see Table [2]). Lower 
panels: Velocity residuals to the adopted fit for 
the two components. 



The "Aa" and "Ab" components correspond to the 
broad- and narrow-line spectra, respectively. 

First, we used clean (Roberts et al. 1987), 
under an adapted version of the original code par- 
ticularly useful for unequally spaced data, for de- 
riving a preliminary orbital period of 147±6 days. 
Next, we applied the Lehmann-Filhes method im- 
plemented in SBOP (Etzel 2004) for a detailed or- 
bital analysis. We computed the orbital solutions 
for periods between Psb2 = 141.0 and 147.0 d in 
steps of 0.5 d and found that the x 2 of the radial 
velocity curves of both broad and narrow compo- 
nents minimized at PsB2 = 143.5 d. We thus as- 
sumed an error of 0.5 d in the orbital period of the 
double-line spectroscopic binary (SB2). The final 
orbital solution is shown in Table [2] and illustrated 
in Fig. H 

The derived spectroscopic orbital period, of 
about 0.39 yr, is roughly 400 times smaller than 
the astrometric orbital period of ^157 yr. Then we 
face a hierarchical triple system containing a close 
SB2 (the components Aa and Ab) with a period of 



Parameter 


Value 


Unit 


•PsB2 


143.5 ± 0.5 


d 


T 


2455452.03 ± 0.18 


d 


e 


0.7834 ± 0.012 




7 


+29.8 ± 0.3 


km s 


M Asi /M Ah 


1.23 ± 0.07 







Aa 


Ab 




K 


75 ± 4 


93 ± 5 


kms 1 


it) 


20.8 ± 0.9 


200.8 ± 0.9 


deg 


a sin % 


132 ± 7 


164 ± 8 


Rq 


M sin 3 % 


9.3 ± 0.8 


7.6 ± 0.7 


Mq 



a few months and a fainter astrometric compan- 
ion (the component B) orbiting a common center 
of gravity with a period of over a century. This 
scenario is also compatible with the good agree- 
ment found between the center-of-mass velocity 
obtained for a Ori Aa and Ab and the systemic 
velocity determined for the cluster by Sacco et al. 
(2008). Following predictions by Hartkopf et al. 
(1996) both quantities should not differ by more 
than ~ 4 kms . 

3. Discussion 

A small region of the spectrum of a Ori Aa+Ab+B 
obtained on 2010/09/09 (corresponding to the 
largest separation between the broad and nar- 
row components) is displayed in Fig. [3J The se- 
lected region contains the He II A4541 A line and 
the Si IIIAA4552,4567,4574A triplet, traditionally 
considered to establish the spectral type of late 
O- and early B-type dwarfs. The spectra of the 
standard stars HD 34078 (AE Aur) and HD 36960 
from the IACOB spectroscopic database (Simon- 
Diaz et al. 2011) are also shown in the figure. 
The two spectra were convolved to the projected 
rotational velocities of the broad and narrow com- 
ponents of crOri Aa and Ab (Section [2]), shifted 
to the radial velocities indicated in Table [TJ and 
diluted by factors 0.500 and 0.275, respectively. 

Three conclusions can be made from Fig.O (i) 
the two spectroscopic components are nicely fitted 
by the spectra of the 09.5 V and B0. 5 V standards, 
(ii) in a lower resolution spectrum (or even at a 
phase in which both components are closer) the 
spectrum of a Ori Aa+Ab+B could be erroneously 
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classi fied as a "single" 09.7V star (e.g., Sota el al. 
20111 ). and (Hi) there is no clear evidence of the B 
component, which should be located close to the 
systemic velocity of the a Ori cluster, making its 
spectral classification difficult. From the observed 
magnitude differences with respect to Aa+Ab (see 
e.g. ten Brummelaar et al. 2000, Caballero 2008a, 
Maiz Apellaniz 2010), the third star must be an 
early-B dwarf, slighly cooler and less massive. 

In spite of the complexity of the system, we can 
estimate approximate individual masses of the hi- 
erarchical triple system by applying the third Ke- 
pler's law (i.e., (M Aa + M Ab + M B )P 2 = a 3 , where 
a = ad is the physical semimajor axis in AU, a is 
the angular semimajor axis in arcsec, d is the he- 
liocentric distance in pc, P is the astrometric or- 
bital perior in years, and Mi are the stellar masses 
in M Q ), the ratio of masses M Aa /M Ah (Tabled]), 
and apparent magnitudes of resolved Aa+Ab and 
B (see above). Using a and P from Hartkopf et al. 
(1996) and assuming d = 385±15pc, we estimate 
the individual masses at M Aa — 19±2, M A \, = 
15.4±1.6, M B = 9±4M , so the triple system con- 
tains over 40 Mq. The system total increases to 
over 60 M Q if we also take into account a Ori C, 
D, and E. 

The existence of an additional powerful ultra- 
violet source and of a deeper central gravity well 
than previously thought has a critical impact on 
the way the astronomers see the a Orionis cluster 
and its surroundings. In addition, the dynamical 
properties of the central, young triple system de- 
serve further observational and theoretical studies 
in the framework of the evolution of highly eccen- 
tric massive binaries and massive star formation 
theories. 

We are in debt to D. M. Peterson for pro- 
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star system. We thank to the anonymous ref- 
eree for a very constructive report, and to J. Maiz 
Apellaniz, R. Barba, and N. Walborn for a care- 
fully reading of the manuscript before publica- 
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Spanish Ministerio de Ciencia e Innovation un- 
der projects AYA2008-06166-C03-01, AyA2008- 
06423-C03-03, and AyA2008-00695. This work 
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[right], and B [center] (measured in the case of the former ones, and expected in the case of the latter one) 
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